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We discuss charmonium production in proton-nucleus collisions at RHIC energies under the assumption oi xf 
and X2 scaling. We find that all the ambiguities due to energy loss are gone at this energy and therefore data will 
reveal the scaling law, if any. These p-A data will also be crucial to interpret nucleus-nucleus data with respect 
to a possible formation of a quark gluon plasma because the extrapolations for charmonium production from the 
present p-A data to RHIC energies, based on the two scaling laws, differ by a factor of four. 



I. SCALING IN CHARMONIUM NUCLEAR 
PRODUCTION 

Charmonium production off nuclei has shown to be 
strongly affected by nuclear effects in proton induced re- 
actions at several beam energies jl]-B]. It is therefore 
natural to wonder whether the nuclear dependence of 
charmonium production is actually driven by a scaling 
variable, independently of the center-of-mass energy y/s 
of the collision. Indeed, a first glimpse at proton- nucleus 
data at SPS |l| and Fermilab tends to indicate that 
J I tp suppression scales with the longitudinal momentum 
fraction Feynman-cc (xp)- 

Such an xp scaling has been predicted long ago within 
the intrinsic charm model ||. Indeed, charmonium pro- 
duction at large xf may be dominated by the soft scatter- 
ing of the intrinsic cc component of the projectile on the 
front surface of the nucleus leading to an unusual A 2 / 3 
dependence at high xf, i-e., a suppression as compared 
to the A 1 dependence of the hard parton fusion process. 
However, the intrinsic charm content of the proton re- 
quired to describe quantitatively the data is somewhat 
larger than what is found from a leading order analysis 
of the EMC charm structure functions ||. 

Furthermore, it is expected from the QCD factoriza- 
tion theorem that nuclear effects at high energy scale 
with the momentum fraction X2 carried by the target 
parton from the nucleus. Furthermore, absorption of 
charmonia in nuclear matter depends on the space-time 
evolution of the expanding cc state as well as on the 
energy dependence of the charmonium hadron inelastic 
cross sections || . Since both the Lorentz factor 7 of the 
cc system in the nucleus rest frame (i.e., which governs 
the expansion of the cc precursor in the nucleus) and 
the incident energy ^Js^h a- re related to the momentum 
fraction X2, final state interactions naturally lead to an 
X2 scaling in charmonium suppression, assuming a given 



mass for the cc state. Similarly, the effects of nuclear 
shadowing Q in charmonium production are expected to 
depend on #2, independently of the center-of-mass en- 
ergy of the collision. However, we would like to point out 
that possible effects of parton energy loss in nuclear mat- 
ter might distort a possible X2 scaling law at low beam 
energy @Jll|. 

On the experimental side, we emphasized that the 
comparison of NA3 § {*/s = 19.4 GeV) and E772 @ 
(a/s = 38.8 GeV) data seems to show that J/ip suppres- 
sion in proton-nucleus collisions scales with the longitu- 
dinal momentum fraction Feynman-x (xf) [§■ However, 
xp scaling is no longer reported in pion- nucleus NA3 
data [@ from ^ = 16.7 GeV to ^/s = 22.3 GeV |J]. The 
X2 scaling assumption appears to be ruled out at small 
x 2 ~ 0.07 when considering NA3 and E866/NuSea p-A 
data, even though a good agreement is found at larger 
X2- Thus, it is our understanding that two quite dif- 
ferent scenarios can be invoked in order to explain the 
present p-A data : i) A universal xf scaling law, valid 
from y/s « 20 GeV on; ii) A X2 scaling law visible at 
higher energy, but distorded in the low energy range by 
effects such as energy loss. Unfortunately, no clear and 
definitive conclusion can be drawn from the present ex- 
perimental results, first of all because considerable error 
bars weaken the confidence one could have in some of the 
most relevant data points, but also because energy loss 
effect is in fact not perfectly understood at NA3 energies. 

In order to distinguish clearly which scaling law, if any, 
governs nuclear effects seen in charmonium production, 
one thus has to turn to higher -^/s, such as presently avail- 
able at the RHIC accelerator. It is the main goal of this 
article to show that energy loss effects are indeed negli- 
gible in this energy range, while X2 and xp scaling law 
assumptions lead to quite different predictions for the 
production rate of charmonium in proton-nucleus colli- 
sions planned in year III of RHIC. Consequences for the 
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interpretation of nucleus-nucleus experiment are eventu- 
ally investigated. 



II. RHIC PREDICTIONS 

Using the 800 GeV proton beam at Fermilab, the fixed 
target experiment E866/NuSea measured J ftp and ip' 
production off several nuclear targets (Be, Fe, and W) 
over a wide kinematical range . The production ratio 



R(W/Be,x F ) = 



9dcr(>W ip(x F ) X) 
184dcr(pBe ij}{x F )Xy 



has been extracted for longitudinal momentum fractions 
— O.f < xp < 0.93. It is found to decrease strongly from 
i?(W/Be)=0.86 at x F w 0.2 down to i?(W/Be)=0.3 at 
x F = 0.9. 

In Q , the x p dependence of the differential production 
ratio has shown to be consistent with 



cx(xf) — 1 



where a{x F ) follows the fitting law 

a E866 (x F ) = 0.960 (1- 0.0519 xp- 0.338 x%) 

in the whole experimental x F range. 

x F scaling. Assuming x F to be the proper scaling vari- 
able, i.e., R independent of y/s, the expected J ftp pro- 
duction rate at RHIC is straightforward. Provided the 
power law <j{A) oc A a is valid even for light systems, 
the ratio R(Au/p,x F ) in p(100 GeV)+A(100 AGeV) is 
simply given by 

R(Au/p,x F ) = 197 aE866(x ^- 1 . 
The simple relation between x F and the J /ip rapidity 
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J ^ = Arcsinh 



Xpyfs 



2 m 



J/4 



determines the rapidity range over which reliable RHIC 
predictions, based on E866 data, can be made. This 
range is summarized in Table |. The corresponding pro- 
duction ratio is displayed in Figure [I] (solid) as a function 
of Ucm '■ Provided Xp scaling, charmonium suppression is 
expected to be small (R ~ 0.8) and approximately con- 
stant in the range — 2 < ytm < 2. The large nuclear ef- 
fects seen in E866 data are found to be relevant at RHIC 
for rapidities greater than j/ C m ~ 

2.6 (cf. Table Q, right 

and Figure |lj). 

X2 Scaling. The scaling with the momentum fraction 
X2 leads to a significantly different suppression pattern. 
Expressing X2 as a function of x F 



Jli> 



4m j/4 



7?) 





full range 


steep drop 


if in E866 


-O.f <x F < 0.93 


0.2 <x F < 0.93 


vH* at RHIC 


-1.9 < yti^ < 4.1 


2.6 < yiL* < 4.1 



TABLE I. x F range covered by the E866/NuSea data and 
the corresponding J/tp rapidity range at RHIC (y/s = 200 
AGeV) (left). Ranges for which the ratio R(W/Be) sharply 
drops are also given (right). 



and the J / ip rapidity as a function of X2 



£!$ = In 



m 



J/4> 



the sudden drop of the production ratio R occurs around 
midrapidity at RHIC energies, as shown in Table [0] 
(right). Furthermore, we note that X2 scaling predic- 
tions, based on E866 data, are limited to the rapidity 
range -1.1 < yllf < 0.8 (Table g, left). 





fall range 


steep drop 


X2 in E866 


0.044 >x 2 > 0.006 


0.028 > xi > 0.006 


yt^ at RHIC 


-1.1 < w^l* < 0.8 


-0.6 < y%+ < 0.8 



TABLE II. x-i range covered by the E866/NuSea data and 
the corresponding J/ip rapidity range at RHIC (yfs = 200 
AGeV) (left). Ranges for which the ratio R(W/Be) sharply 
drops are also given (right). 

Thus, the steep fall of charmonium production seen 
in E866 data appears in two distinct rapidity ranges 
at RHIC, depending on whether it scales with if or 
X2 ■ Consequently, the rapidity dependence of charmo- 
nium production measured in proton-nucleus collisions 
at RHIC may allow one to disentangle these two assump- 
tions. 

The Model. The expected charmonium rate at RHIC 
under the X2 scaling assumption is discussed in this sec- 
tion. In order to extrapolate to x 2 ranges not covered by 
E866 data (Table ||) , we use a recently advanced X2 scal- 
ing model based on absorption mechanisms [ p^[ . Here, it 
will be complemented by a calculation of parton energy 
loss, that has been proven to be relevant at low beam 
energy @. 

The model is based on a standard two stage color neu- 
tralization scenario. A compact partonic (cc)$ state is 
assumed to be produced by parton fusion in a color octet 
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state. The latter will then turn into a singlet (cc)i state 
after a color neutralization time rg_»i. During the first 
stage, the octet state may interact with surrounding nu- 
cleons with an octet cross section a^ c ^ sN , taken to be 
independent of its transverse size r cS . Adopting an en- 
ergy dependence related to J/ip photoproduction data, 
the octet nucleon cross section reads as 



&{cc) s N = °~8 ' 



lOGeV 



0.4 



The singlet-nucleon total cross section cr^ cS ^ lN is propor- 
tional to the square of the radius r cS and has a same 
energy dependence. We assume that the transverse dis- 
tance between the quarks increases linearly with time 
until the respective hadronic radius ri is reached. In the 
cc rest frame, we write 



ro + v c£ t 

n 



if r C c(r) < n, 
otherwise, 



where i stands for J/ip, ip' and \- Therefore, the time 
dependent cross sections for singlet states are given by 



r (ce) 



lOGeV 



/ r- c5 (r ) x 2 



Using one common set of parameters fitted to 
E866/NuSea results]], almost all other available data on 
charmonium production in p-A collisions can be repro- 
duced within this approach Jl3| . 

As this inelastic cross section is intimately related to 
the time evolution of the cc system that travels through 
the nucleus, charmonium absorption within this model 
is driven by the Lorentz boost 7 of the cc pair in the 
nucleus rest frame: 



7 



rn c 



2x2m p 



where xi is the momentum fraction carried by the target 
parton and m cS ~ mj/^. It therefore leads to X2 scaling 
for the nuclear charmonium absorption. 

Let us now turn to the predictions of the model. Us- 
ing the same set of parameters, we have calculated char- 
monium suppression in p(100 GeV)+A(100 AGeV) colli- 
sions over a wide ra-nge of rapidities ( — 3 ycm 

< 3). The 

J/ip (dashed) and ip' (dotted) production ratios i?(Au/p) 
are depicted in Figure [l]. As expected, both ratios dra- 
matically drop around mid-rapidity. In our model, this 
fall is a consequence of the increase of color-octet nucleon 
inelastic interactions (as compared to the color-singlet 
nucleon interactions) with 7, i.e. with xf- Nevertheless, 
we would like to remind the reader that such a qualitative 



r = 0.15 fm, « C c= 
a s = 22.3 mb. 



1.85, 0-1 = 2.1 mb, r 8 ^i= 0.02 fm, 
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FIG. 1. Ratio R(Au/p) versus y cm of produced J/i/)'s 
within the E866/NuSea fit (solid line) in p(100 GeV) + 
Au(100 GeV) collisions, assuming if scaling. The model pre- 
diction for J/tp (dashed) and ip' (dotted) production ratio is 
also shown. The arrows show the dimuon arm acceptance in 
the PHENIX experiment. 



feature can also be seen as a direct consequence of the x^ 
scaling assumption, together with the particular xi de- 
pendence of the E866/NuSea data. Actual predictions of 
the model are in fact the behavior for the rapidity ranges 
Vcm < -1 and y cm > 1. 

At large y CI[L (i.e., large 7 factor), the partonic cc state 
travels through the entire nucleus in a color octet state 
and hence leads to a constant i?(Au/p)Q as low as R w 
0.2. 

The smaller the values of y cm the more likely the color 
neutralization to happen inside the nucleus. For neg- 
ative y cm the charmonia have time to expand to their 
asymptotic size inside the nucleus. Then, J/ip and ^'s 
have different cross sections with the surrounding matter 
and consequently the respective production ratios differ 
(Fig. Q). Finally, let us also mention that these general 
trends are found to be similar in lighter systems. 



III. SHADOWING AND ENERGY LOSS 

So far, charmonium suppression has been attributed to 
the inelastic interactions of the cc pair in nuclear matter. 
Let us investigate in this section the effects of both shad- 
owing and parton energy loss in the nuclear dependence 
of charmonium suppression at RHIC energies. 



2 Actually a slight decrease is expected due to the weak en- 
ergy dependence of the octet-nucleon cross section (see sec- 
tion II.A) 
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FIG. 2. Ratio R(Au/p) versus j/ cm of produced J ftp's 
in p(100 GeV)+ Au(100 AGeV) collisions predicted by the 
model assuming nuclear absorption only (solid), nuclear ab- 
sorption with EKS98 shadowing (dashed), nuclear absorption 
with BDMPS energy loss (dotted). 



Shadowing. It is known from the first EMC data that 
structure functions (x) in a heavy nucleus strongly dif- 
fer from that measured in light targets. At small x-i <C 1 
(i.e., large yiia), the depletion of parton distributions 
(shadowing) lowers the charmonium nuclear production, 
hence leading to a production ratio smaller than one. 
On the contrary, possible gluon anti-shadowing between 
0.1 < xi < 0.3 enhances charmonium production in a nu- 
cleus in the range —3 < yim < —1.8 at y 7 ! = 200 GeV. 

Shadowing effects on J/ip suppression in proton- 
nucleus collisions have been investigated by several au- 
thors H . They all report that it leads to a J/ip suppres- 
sion at large xp at Fermilab energy in agreement with the 
general trend reported by E772 {§ and E866/NuSea (§. 

We compute in Figure ^ (dashed) the J/ip produc- 
tion ratio in p(100 GeV)-Au(100 AGeV) assuming both 
nuclear absorption and shadowing effects. The produc- 
tion cross sections are given by the color evaporation 
model with MRST parton distributions 14 1. Nuclear 
shadowing is evaluated using the parameterization of 
deep inelastic scattering and Drell-Yan data by Eskola 
et al. (EKS98) 0. A glimpse at Figure | indicates 
that the general trend predicted for the rapidity de- 
pendence of J/ip suppression still holds when shadow- 
ing of target partons is taken into account. We also 
note that measurements in the negative rapidity range 
(—2.2 < yim < —1.1) may signal gluon antishadowing 
effects. Furthermore, charmonium production at large 
rapidities (ytl/t < 2) will reveal a possible saturation of 
shadowing at small X2 ~ 10~ 3 . 

While the EKS98 shadowing based on the factoriza- 
tion of the nuclear structure functions leads only to a 
small suppression as compared to nuclear absorption, 
it is worth pointing out that a rapidity dependence for 



J/ip suppression at RHIC similar to the one in Figure || 
has also been predicted as coming from coherence effects 
only @. 

Energy loss. The incoming parton from the projectile 
may suffer multiple scatterings while traveling through 
the nucleus A, hence leading to a shift of its momentum 
fraction from x\ + Axi to x\ at the time of the cc produc- 
tion. In addition to final state interactions, ip production 
in a nucleus will thus be suppressed roughly by a factor 

cr(pp -» iP( Xl + Axi) X) fa{pp -* ip(xi)X) (1) 

Because of the dramatic drop of the parton distributions 
at large x (especially that of the gluon) , one may reason- 
ably expect a significant suppression even for small shifts 
of the momentum fraction x\ . 

A lot of theoretical effort has been devoted in the last 
few years to relate the momentum fraction shift Ax\ to 
the atomic mass number A and the center-of-mass energy 
\fs of the collision |p7|-|l9|. For the numerical computa- 
tions to follow, we shall adopt the model by Baier et al. 
(BDMPS) jl9| which assumes the radiative energy loss 
per unit length AE/ Az to depend linearly on the length 
z oc A 1 / 3 covered by the parton in the nucleus, i.e., 

AE Al/3 

— — oc A 1 ' 3 
Az 



which leads to a momentum fraction shift 
Axr = -A 2 / 3 



(2) 



where k is a free parameter^. While the momentum frac- 
tion shift (|^) vanishes at large energy y/s, energy loss ef- 
fects on charmonium production are expected to be size- 
able at low incident energies when Ax± becomes large. 
We already emphasized that such a mechanism may thus 
be responsible for the lack of x^ scaling seen from NA3 
to Fermilab data in J/ip production [pl}|Jll[ . 

Let us now estimate more quantitatively the effects 
of energy loss in J/ip production at RHIC energy. We 
apply the BDMPS model jll| with a quark energy loss 
AE/Az = 0.3 GeV/fm (k = 0.086 GeV 2 ) in a gold 
nucleus^. Energy loss of this magnitude together with 
the absorption model appears to be consistent with both 
NA3 and E866/NuSea J/ip data @. Again, the pro- 
duction cross sections in (Q) are computed within the 
color evaporation model and using MRST parton distri- 
butions Q. 

The J/ip production ratio i?(Au/p) in p(100 
GeV)+A(100 AGeV) collisions including energy loss is 



We assume that the perturbative analysis of BDMPS is 
justified at the energies of RHIC. This may no longer be true 
at SPS energy where the incoming parton energy is small. 

4 The corresponding gluon energy loss is expected to be 
dE/dz\ g = 9/AdE/dz\ q « 0.7 GeV/fm in a gold nucleus. 
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displayed in Figure || (dotted). Energy loss proves to 
be relevant in the nucleus fragmentation region and to 
be negligible otherwise. Therefore, the large suppression 
predicted at RHIC by the 22 scaling assumption is not 
altered by such a mechanism. We further note that pos- 
sible energy loss effects may be seen and better under- 
stood [[y]] at HERA-B where charmonium and Drell-Yan 
production will be measured in p(920 GeV)-nucleus col- 
lisions for -0.5 < x F < 0.3 (20). 

Whereas we focus here on hadron-nucleus reactions, we 
would like to emphasize that parton energy loss in hot 
QCD matter (or "jet quenching", |21|) is expected to be 
even much larger than that in nuclei of about one order 
of magnitude We note that the PHENIX collab- 

oration at RHIC recently reported on the suppression of 
the yield of high px hadrons that may be attributed to 
finite parton energy loss |2^| . In view of that, determi- 
nation of finite energy loss in proton-nucleus collisions 
would certainly help to constrain estimates obtained in 
nucleus-nucleus collisions. Finally, let us also mention 
that elastic scatterings of the produced cc pair with the 
nucleus (final state energy loss) may also lead to large 
effects in the nuclear dependence of J/ip production p3[ , 
although the radiative energy loss of heavy quarks may 
be much smaller than the one of light quarks (24] . 

Furthermore, we have checked that this particular ra- 
pidity dependence is expected to persist when nuclear 
shadowing of target partons, predicted by the EKS98 pa- 
rameterization Eql, is taken into account. 



range may allow to constrain formation times estimates 
as well as the time evolution of the partonic system. Both 
are not precisely known so far. 

One may wonder whether energy asymmetric collisions 
like p(250 GeV)+A(100 AGeV) j|| allow to investigate 
scaling with the present detector acceptance. As shown 
in Figure I this is not very promising. In order to observe 
the steep drop of R in the acceptance region, one has in 
fact to lower the gold beam energy even further, e.g., 
p(250 GeV)+Au(65 AGeV). 
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FIG. 3. Idem as Figure 1 in asymmetric p(250 GeV) + 
Au(100 AGeV) and p(250 GeV)+ Au(65 AGeV) collisions. 



IV. EXPERIMENTAL CONSIDERATIONS 

The results presented so far show that charmonium 
production in p-A collisions at RHIC energies will allow 
to discard experimentally one of the two presently ad- 
vanced scaling assumptions in charmonium production. 
For ruling out one of the assumptions it is indeed suffi- 
cient to measure charmonium production in the rapid- 
ity interval 1 < y cm < 2. Such a measurement can 
be carried out by the present setup of the PHENIX 
experiment which is sensitive to charmonium produc- 
tion in the dimuon decay channel for the rapidity range 
1.2 < |yiab I < 2.4, represented by black arrows in Fig- 
ure [|. 

However, if one would like to confirm scaling, either 
in xf or X2, a precise measurement in the region where 
R drops sharply is necessary, i.e. either in the rapidity 
range 2.6 < yH^ < 4.1 or in —0.6 < y'lm < 0.8 respec- 
tively (see Tables | and |l[ right). These domains are 
most unfortunately not covered by the dimuon spectrom- 
eter acceptance, although the midrapidity region could 
be explored via the e + e~ decay channel of J/ip. 

As for the backward dimuon spectrometer, it will per- 
mit to investigate formation time effects in charmonium 
production, as demonstrated in Figure [l]. In particular, 
the ip' over J/ip ratio in various nuclei in this rapidity 



V. BASELINE TO NUCLEUS-NUCLEUS DATA 

Heavy quarkonium suppression is expected to signal a 
quark-gluon plasma (QGP) formation in nucleus- nucleus 
collisions p^ j. On the contrary, no deconfinement is ex- 
pected in proton- nucleus collisions. Therefore, p-A mea- 
surements provide a natural benchmark to which a pos- 
sible additional suppression due to QGP formation has 
to be compared with. 

However, we have seen that the expected J/ip rate 
in proton-nucleus collisions is strongly dependent on the 
scaling law of the underlying mechanism responsible for 
charmonium suppression. To illustrate this, let us give 
an estimate for J/ip production ratio in inclusive Au- 
Au collisions at y/s = 200 AGeV, assuming no abnormal 
suppression (i.e. no QGP formation). In this case, the 
production ratio at a given y cul is expected to be roughly 

R(Au - Au, y cm ) ps R(Au - p, y cm ) R(Au - p, -y C m) 

Assuming xp scaling, we expect a constant production 
ratio R(Au/p) ps 0.64 in y cm e [-2,2]. Provided x 2 
scaling prevails in charmonium nuclear production, this 
ratio is also constant and as low as R(Au/p) pa 0.16. 
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Thus, the extrapolation from presently available p-A 
data to the hadronic suppression expected at RHIC en- 
ergies yields an uncertainty of a factor of four. There- 
fore, an unambiguous judgement whether J/ip suppres- 
sion seen in A-A collisions at RHIC is due to QGP for- 
mation has to await for p-A measurements at the same 
energy over a wide range of rapidities. 
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